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ABSTRACT 

High-contrast imaging instruments are now being equipped with integral field spectrographs (IFS) to 
facilitate the detection and characterization of faint substellar companions. Algorithms currently en- 
visioned to handle IFS data, such as the Locally Optimized Combination of Images (LOCI) algorithm, 
rely upon aggressive point-spread-function (PSF) subtraction, which is ideal for initially identifying 
companions but results in significantly biased photometry and spectroscopy due to unwanted mixing 
with residual starlight. This spectro-photometric issue is further complicated by the fact that algo- 
rithmic color response is a function of the companion's spectrum, making it difficult to calibrate the 
effects of the reduction without using iterations involving a series of injected synthetic companions. 
In this paper, we introduce a new PSF calibration method, which we call "damped LOCI", that seeks 
to alleviate these concerns. By modifying the cost function that determines the weighting coefficients 
used to construct PSF reference images, and also forcing those coefficients to be positive, it is possible 
to extract companion spectra with a precision that is set by calibration of the instrument response 
and transmission of the atmosphere, and not by post-processing. We demonstrate the utility of this 
approach using on-sky data obtained with the Project 1640 IFS at Palomar. Damped-LOCI does not 
require any iterations on the underlying spectral type of the companion, nor does it rely upon priors 
involving the chromatic and statistical properties of speckles. It is a general technique that can readily 
be applied to other current and planned instruments that employ IFS's. 
Subject headings: exo-planets, instrumentation, speckle calibration, high contrast imaging 



1. INTRODUCTION 

Current instruments dedicated to the imaging of exo- 
planets have a contrast floor that is set by the imper- 
fections on the surfaces of the the optical train. In the 
case of a space observatory, as show n for instance by 
Krist et al.[ fl2005| ); |Kalas et ah] ( |2008p , when the opti- 
cal surfaces corrugations are stable in time, the instru- 
ment response can be calibrated and subtracted. How- 
ever, in the ground based regime, the PSF varies over 
timescales that r ange from minutes to hours e .g |Hink-| 
ley et al.| f2QQ7[ ); |Fitzgerald fc Grahaml ( [2006J ). Since 
these so called quasi-static speckles are the limiting fac- 
tor of ground and spaced-based observations with re- 
gard to faint companion detect ability, their properties 
have been studied extensively over the past few years. 
It was found that their theoretical properties ( |Perrin 
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l et al ||2003| ;|Aime fc Soummer ( 2 004| ) ; |Soummer fc Fer- 
|rari| (|2007|) ; [Soummer et al. ( 2007a JTmatch the statistics 



observed~bri -sky (Fitzgerald fc Graham (2006); Hinkley 



eta^ ( |2QQ7t ). 

Several observing techniques and data reduction algo- 
rithms have been devised over the past few years in or- 
der to optimize the sensitivity of current imaging instru- 
ments to faint exo-planets. Observations are generally 
carried out by acquiring of a series of images in which the 
observer deliberately chooses to introduce some diversity 
in the imaging method. The purpose of this diversity is 
to force the image location of the residual speckles to vary 
from frame to frame, while the image of the sky is immo- 
bile. This is the case in t he Angular Diff erential Imag- 
ing (|Marois et al.| fl2006D; [Leconte et al. | (|2007 D) , Spec- 
tral 



tral Different ial Imaging (Racine et al. ( |l999p ; |Marois| 
et al.j fl2004D;|Biller et al | OT)7|) ), chromaWHifferential 



imaging (Sparks fc Ford (2002)) and Coherent Imaging 
( |Guyonet al.| ( |2010| )). Post-processing techniques involve 
a manipulation of this series of images whose purpose is 
to remove the artifacts due to speckles and leave intact 
the signal of a companion. The powerful combination 
of ADI data processed using a Locally Optimized Com- 
binatio n of Images (LOCI) PSF subtraction algorithm 
allowed pVIarois et al. (2008) to obtain one of the first im- 
ages of a planetary system . The LOCI PSF subtra ction 



technique was devised by Lafreniere et al. (2007) and 



is quickly becoming a standa rd in exo-plane t searches 



(|Marois et aT| (|2008|); rT halma nn et al.| fl2009| ); |Metchev 



et al.| ( |200^ " |(Jrepp et al.| ( |2010[ )). Other techniques cur 



rently discussed in the literature involve a-priori models 
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of the PSF and the telescope Mugnier et al. (2009 ); Burke 
fc Devaney (2010) or spectral deconvolution based on 
priors about the u nderlying spectrum of the companion 
Mesa et al. ( 2011[ ). On the contrary, LOCI only requires 
minimal priors about the behavior of the instrument and 
builds a composite reference PSF solely on the data ob- 
tained during observations. 

In a previous paper we reported a recent implementa- 
tion of LOCI that is applied to process data from an Inte 
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gral Field Spectrograph (IFS), Crepp et al | ( |2010| ). Our 
results were obtained using Project 1640. The lull system 
is composed of the Palomar Adaptive Optics (PA L AO) 
followed by a stellar coronagraph ( jSoummer (2005)), and 
the science science IFS. Further d e tails a bout the mstru- 
ment are found in Hinkley et al.| (2011). The raw data 
from the infrared detector is calibrated and re- arranged 
in order to produce a data cube of images at 23 differ- 
ent wavelengths across the J and H bands (R ~ 40). 
The spatia l and spectral data extract ion procedures are 
detailed in Zimmerman et al. ( 2010a[ ). The P1640 con- 
figuration prevents the held rotation that is necessary to 
carry out ADI observations, but allows both chromatic 
and temporal diversity. Using an IFS behind a high con- 
trast starlight suppression system is key to the two main 
science objectives of P1640: (a) detection of faint com- 
panions around nearby stars, (b) in-situ low resolution 
spectral characterization of any companions discovered. 

The detection of faint point sources is facilitated by the 
chromatic diversity of the IFS and allows the observer to 
apply advanced static speckles calibration methods. In- 
deed, if all the errors are in an instrument pupil plane, 
diffractive optics predicts that the speckle pattern will 
stretch with wavelength, while the scale of the astronom- 
ical image is independent of wavelength. In an IFS, the 
series of measurements of the quasi-static speckle pat- 
tern is simultaneous. This alleviates the problem of cor- 
rugations evolving temporally between sequential obser- 
vations. On the P1640 camera the speckles located at 
1" are at the 5 x 10 -3 level and de-correlate in inten- 
sity by roughly 1 percent per minute. This implies that 
purely sequential observations, separated by 5 minutes 
in time, can at best reach a contrast of 2.5 x 10 -4 after 
subtraction. Using an IFS alleviates this problem as it 
provides reference PSFs which are contemporaneous re- 
alizations of the speckle pattern present in the science 
exposure. Each wavelength slice can be interpolated in 
order to produce rescaled data cubes in which the speckle 
pattern's scale remains constant. In this rescaled space 
the companion moves radially. If the variations of the 
speckle's brightness as a function of wavelength is know 
a-priori, then these rescaled images can be subtracted 
one to another to detect a faint companion . Thi s ap- 
proach was first discussed in |Sparks fc For d ( 2002 ) and 
has been suc cessfully applied to both s imulated and on- 



Thatte et al. 



(2007) . It 



sky IFS data |Vigan et aL| ( |2QQ8| >;f 
is very powerful as it uses reference images that have all 
been acquired simultaneously, thus reducing the impact 
of the temporal evolution of the speckles to time sc ales 
shorter than the exposure time. In lCrepp et al.| ( |20To| ) we 
presented an algorith m that combines th e chromatic dif- 
ferential imaging in |Sparks fc F ord (2002) with the LOCI 
PSF sub traction approach. Wh ile LOCI was first pre- 
sented in Laf reniere et al.| ( |2007[ ) as a method to reduced 
ADI images, P1640 does not take advantage of field ro- 



ta tion. Instead , the b aseline P1640 reduction presented 
in Crepp et al. (2010) relies on the IFS' radial chromatic 
stretch ol the PSF and on temporal diversity. It inte- 
grates them both in a LOCI-type reduction scheme. We 
showed that this quasi-static speckles calibration method 
provided a gain greater than an order-of-magnitude in 
the contrast between raw and processed images, and an 
ultimate detection level ~ lx 10 -5 at 1" in a 20 minutes 
H-band exposure on a three magnitude star. Further 
improvements in contrast, while maintaining reasonable 
exposure times, require cohere nt wavefront calib r ations 
such as the ones presented in Borde fc Traub| ([2006 



Guy on et al | fl2010| ) ; |Give 7 on et al.| ( |2007p ; (Pueyo et 
( 2009 ]), lor the ca s e of space- based observ ations and Wal- 



36); 
. al. 



|lace et aL] ( 2009 ) ; Sauvage et al. (2007) for the case ol 
ground- based telescopes. 

A second advantage of the P1640 IFS is that, once a 
companion is detected by means of the aforementioned 
chromatic diversity, its spectral information is an imme- 
diate by-product. Over the past ten years, IFS have been 
used to ch aracterize companions with favorable contra st 
ratio (e.g. |McElwain etaL] ( |2007[ ); [Patience et aT](|2010| )). 
When the companions is fainter than the speckl e floor, 
then reductio ns such as the ones presented in |Crepp| 



et al. (2010) ought to be carried out to first detect 
the companion. However, blind applications of LOCI 
post-processing tailored for detections to spectral data- 
cubes present the major caveat of altering the spectro- 
photometry and possibly introduce artificial features in 
the low resolution of the detected object. Recent results 
obtained using the OSIRIS IFS at the Keck observatory, 
show a depletion of sixty percent of the co mpanion flux 



that is constant over the spectral window (Bowler et al. 
k). In this example the depletion is gray and as 



such has little impact on the scientific conclusions of the 
observations. However PI 640 seeks for faint objects at 
close angular separations whose signal is buried within 
the bright speckle field, thus requiring more aggressive 
PSF subtraction. Moreover PI 640 entirely relies on chro- 
matic diversity, which can potentially lead to a poor con- 
ditioning of the inverse problem associated with LOCI 
(see below). In practice, when testing our reduction on 
synthetic objects injected within the control radius of 
the AO system, we observe that the post-LOCI retrieved 
spectral information is considerably biased with a non- 
gray algorithmic response. This biasing of the spectral 
information is entirely due to the PSF subtraction algo- 
rithm, and can be an order-of-magnitude larger than the 
next most important uncertainty, due to the telluric cali- 
bration described in Zimmerman et al. (2010a). It stems 
from two independent phenomena: a bias caused by par- 
tial fitting of the companion flux with the starlight in the 
reference frames, and a spectral cross-talk between actual 
companion flux between neighboring channels. The com- 
bination of these errors is very often the dominant source 
of uncertainties in spectroscopic measurements carried 
out using post-LOCI subtracted IFS data-cubes. 

Herein we introduce a new method that lowers the er- 
ror in the measured Spectral Energy Distribution due 
to the reduction algorithm close to the level of uncer- 
tainty in the telluric calibration. We first introduce a 
modified version of the LOCI algorithm that alleviates 
the problem of flux depletion of faint companions that is 
common to most locally optimized PSF subtraction al- 
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gorithms. The "damped LOCI" algorithm introduced in 
this communication seeks to maximize the residual flux 
of a putative companion while minimizing the noise due 
to the quasi-static speckles, whereas standard LOCI im- 
plementations only focus on the latter objective. This 
is achieved by re-formulating the least-squares problem 
associated with finding a locally optimal composite refer- 
ence PSF: we augment the cost function associated with 
the determination of the composite PSF coefficients with 
a penalty term that scales with the flux of a potential 
companion. "Damped LOCI" thus solves the following 
optimization problem "minimize the least-squares fit to 
the local speckle field while maximizing the least-squares 
fit in the area where a known companion is presenf . 

The formalism and the notations underlying LOCI and 
"damped LOCI" are details in appendix A. §2 focuses on 
the specifics and the implementation of each algorithm. 
Both approaches yield comparable performances in terms 
of band-averaged faint companions' SNR. Appendix B 
discusses the second order impact of each method on the 
detectability problem, and in particular the advantages 
of "damped LOCI" when the inverse problem associated 
with the PSF reconstruction is ill-conditionned. In §3 we 
discuss the core matter of this paper, namely the mea- 
surement of the spectrum of faint companions. When 
companions have already been detected we show that 
a series of binary PSF weightings, or masking, can be 
used in conjunction with classical LOCI methods to re- 
trieve the spectrum of relatively bright objects. We il- 
lustrate however that the performance of such methods 
degrades rapidly with increasing contrast, and show that 
using masking in conjunction with "damped LOCI" de- 
couples spectral precision and contrast. We show that 
this method yields in all cases better absolute and rel- 
ative spectro-photmetric accuracy when compared with 
classical LOCI reductions. In particular, for objects as 
bright as the residual speckles, we obtain estimates of the 
Spectral Energy Distribution of the order of the telluric 
calibration uncertainties of the spectrograph . For faint 
objects "damped LOCI" yields SED uncertainties consis- 
tent with higher SNR (less "photon starved " ) state of the 
art band-averaged photometr i c estimat es (Marois et al 



(2008); 



Thalmann et al. 



Metchev et al 



(2009p^ 
P1640 



( |2QQ9D ; 

the architec ture ol future exo- planet imagers an d 
are similar (Graham et al. (2007); Dohlen et al. (20061)) 
the results presented here are broadly relevant to the tu 
ture exo-planet imaging endeavors. 



pixels, the subtraction zone ,<S zone. Crepp et al. M2010) 
presented how such an approach can be utilized to re- 
duced IFS data. Appendix A details the specifics of the 
frame selection and the notations underlying the present 
manuscript. In this section we focus on the kernel of the 
LOCI algorithm, namely the inverse problem associated 
with the determination of the optimal reference weights. 

The search for the coefficients of the composite refer- 
ence is based on a least square cost function integrated 
over the zone O. 



mm {c k } 



dxdy w(x,y) 



-i 2 



T(x,y) -^c k R k (x,y) 



£ 

k 



c k Rk\ 



where T stands for the target image, {R k } is the en- 
semble of reference images, {ck} the ensemble of LOCI 
coefficients, w(x,y) is a weighting function that can be 
used to mask certain pixels, Nn the number of reference 
frames in 7£, and ||.|| w ,e> denotes the L 2 norm, weighted 
by w, over the optimization zone. In matrix notation we 
re- write Eq. [T]as: 

||r-5>Jfclko = c r M^ c-2V<?' T > c+S^ 

k 

(2) 

(R R) 

where c is a Nn dimensional vector, M^, ' is an N-jz x 

Nfc matrix, V^f is a N dimensional vector, and S^' T ^ 
a scalar. They are defined by: 



ci,c 2 , ...,Cfc, ...c^ 



(3) 



M ( o ,R) [p,q]= dxdy w(x,y)R p (x,y) R q (x,y) (4) 
Jo 

v cf ' T) [p] = / dxdy w(x, y)Rp(x, y) T(x, y) (5) 
Jo 



5 e> 



o 



dxdy w(x, y)T(x, y) T(x, y) (6) 



Using these notations, the quadratic residual minimiza- 
tion problem associated with standard LOCI implemen- 
tations can be solved by taking the derivative of of 
Eq. |2lvith respect to c, 



2. STANDARD LOCI AND DAMPED LOCI 
IMPLEMENTATIONS 

2.1. Standard LOCI 

The Locally Optimized Combination of Images con- 
sists of reconstructing a reference PSF from IFS data. 
It is based on an inverse problem which seeks minimize 
the residuals between a given astronomical image, the 
target image to be reduced, a nd a calibration image. 
The algorithm, as introduced by |Lafreniere et ah] ([2007 ), 
seeks to find an optimal combination of images, within 
an ensemble of reference frames, that locally minimizes 
the least-squares residual between the target frame and 
the composite reference, in a large optimization zone of 
the image, O zone. This weighted reference is then sub- 
tracted from the reduced image for only a small subset of 
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{Ck} 



(R,R) 
O 
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(7) 



If the ensemble of references, {R^ £ 7^}, forms a com- 

( R R) 

plete basis set of the (9-zone, then the matrix M^, ' is 
invertible. If the speckle diversity is only temporal, and 
if each frame corresponds to an independent realization 
of the speckle pattern, then the inverse problem in Eq. [7] 
is well conditioned when the number of PSF cores in the 
optimization zone matches the number of reference im- 
ages, Na — In this case the optimal weights are 
defined as 



cloci = (Mg^)- 1 ^'^ 



(8) 
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However when using IFS data, the color slices within 
each exposure correspond to the same realization of the 
speckle field, albeit seen at different wavelengths. As a 
consequence they are highly correlated and this results 

(R R) 

in poorer conditioning of M^, ' . This conditioning is a 
measure of how well Na matches the number of linearly 
independent references frames in the zone O. 

Carrying out a selection of the most relevant slices in 
7£, so that the number of independent references matches 
Na, the numbe r of PSF cores, oug ht to be done locally 
and iteratively ( |Marois et al. (2010)). It represents a con- 
siderable computational overhead and in practice, LOCI 
reductions are carried out with Na constant over large 
portions of the image. Then, in the zones for which the 
problem is ill-posed, generic pseudo-inverse routines or 
eigenvalues truncations provide a robust alternative, at 
the cost of losing track of the physics in the process. 
Indeed, such inversion methods lead to the subtraction 
of the least correlated modes in the reference images. 
If the flux of the companion in the target frame can 
be described using some of these independent reference 
modes, then some of its flux is subtracted by the com- 
posite reference. The photometric information is then 
biased. Choosing an optimization zone much larger than 
the subtraction zone, a larger Na, mitigates this issue. 
It increases the contribution of pixels that do not see 
any companion flux and reduces the likelihood of fit- 
ting modes associated with the companion with one of 
the reference modes. However this occurs at expense 
of loosing some of the local characteristics of the sup- 
pression. Such a procedure can be efficient outside of 
the control radius of the AO system, where the speckle 
halo is relatively uniform ( |Bowler et al.| ( |2010| )). However 
when the speckles are bright and structured, optimizing 
the LOCI parameters in order to obtain the best spec- 
tral retrieval quickly becomes an intractable problem. In 
other words, when the minimization is ill-posed, pseudo- 
inverse routines aim at minimizing the speckle noise in 
the reduced image without considerations of the flux of 
a potential companion. Whereas such an approach is 
adequate for faint companion detection, it can be prob- 
lematic when quantifying the spectro-photometry of the 
discovered point sources. We now explore a modification 
of the LOCI that resolves these issues. 

2.2. Flux preserving locally optimized speckle 
suppression: damped LOCI 

"Damped LOCI" relies on the same zonal approach as 
above. The modification we introduce resides in posing 
differently the inverse problem that defines the modal 
coefficients of the composite reference. Instead of solely 
minimizing the residual least-squares fit between target 
frame and composite reference in the optimization zone, 
we constrain the possible combinations of coefficients 
over which the solution may be sought. We first impose 
positivity, c k ^ for all fc, in order to avoid the 
pathological case where two adjacent coefficients are 
large in absolute value but of different signs. Such a 
situation would lead to high spectral cross-talk from 
the companion flux at neighboring wavelengths. While 
this condition might appear too stringent, we show 
below that ensuring the positivity of all the reference 
coefficients allows us to introduce concisely in the 



minimization problem a penalty term that scales with 
opposite of the flux a companion. Because we are 
interested in devising a reduction method that does not 
modify the spectro-photometry of a detected compan- 
ion, the quality of the fit in the zone S is of critical 
importance. Ideally, even if the flux from a companion 
is present in the zone <S, we want the algorithm to 
only fit the stellar contribution in that zone, leaving 
the companion contribution unperturbed. In essence 
"damped LOCI" seeks composite references that are 
the result of a competition between minimizing the fit 
residual in the zone O and maximizing it in the zone S, 
under the constraint that all the composite coefficients 
are positive. Since the zone O is much larger than <S, 
the minimization of its residual fit should prevail, and 
as a consequence the speckles will be subtracted in a 
way that preserves more of the companion flux than in 
classical LOCI implementations. 

Formally, the squared difference between T and the 
composite reference integrated over S can be written as: 



4 2) = | \T - cR\ U s = c T M^ R) c - 2Vf T) c + sf T) 

(9) 

where we have replaced the O subscript by the S sub- 
script for the matrix and vectors. When there is flux 
from a companion in the S zone, we want this quantity 
to remain large while the residual in the zone O is min- 
imized. We can thus use an augmented cost function to 
write the optimization problem as follows: 

min {ckyA} { J (c fe , A)} (10) 
min {ckA} {c T M^ R) c - 2Vg*' T) c + S^ T) (11) 

-A(c T 



^c-2Vf T) c + sf T) )} 



where A is multiplicative constant that al- 
lows us to include the flux preservation con- 
straint in the LOCI minimization. The term 
-A(c T M { f' R) c - 2V^' T) c + SJ' T) ) is called 
the flux preservation penalty term and —A is the 
Lagrange multiplier associated with this term. If one 
seeks to maximize the flux of a potential companion in 
the S zone by minimizing Jo(cfe,A), then A ought to 
be positive. As a consequence, this cost function is not 
positive definite, making the practical implementation 
of this optimization difficult. Forcing positivity for the 
Cfc coefficients leads us to ensure the positivity of the 
second term in Eq. 10 and thus facilitates the search for 
the optimal coefficients. 

To do so, we project the LOCI processed image, 7 P , 
over all the modes present in the collection of references. 



r=T-cR 



(12) 



£ c k I wR k I p = c = Vg^ c - c T M^») 

where we have taken the inner product over the S- zone 
of the composite reference with both sides of Eq. 12 
Replacing the quadratic term in Eq. [9] by the quadratic 
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term in Eq. 13 yields 



if = s 5 



(T,T) 



Yf' T) c. 



(14) 



Since the condition c& ^ for all fc is being enforced 
and the actual images only include positive counts, then 

(T ( R T} 

S s ' > and ' c > 0. In order to devise an aug- 
mented cost function that consists of the sum of positive 

(2) 

scalar quantities, we need to show that maximizing Ig 

( R 

is equivalent to minimizing of ' c. When there is no 

companion in the S zone then JV^' 7 ^ c\ <C IV^'" 7 ^ c|, 
since the flux in the processed image is smaller than the 

flux in the target image. Even if V^' 7 ^ c happens to be 
negative, which means that the subtraction worked too 
well, the integrated square flux in the subtraction zone 

( R 

is dominated by the behavior of ' c. When there 
(R I p ) 

is a companion then ' c > and has the same 

sign as V^' 7 ^ c. Thus the integrated squared residual 
in the subtraction zone behaves as the opposite of the 

positive scalar V^' 7 ^ c. As a consequence we modify 
the augmented cost function shown in Eq. [lOj in order 
to avoid the use of a negative Lagrange multiplier asso- 
ciated with the full quadratic form of the least-square fit 

in the S zone. Instead we choose to use V^'" 7 ^ in con- 
junction with a positive Lagrange multiplier as a penalty 
term. 

Since V^' 7 ^ scales as the opposite of the residual fit in 
the *S zone, this modification re-formulates the problem 
of minimizing Eq. [10] in a tractable way. It allows us 
to write the following positive definite augmented cost 
function: 

J, (c, A) = c T Mg' R) c - 2V(?' T) c + S<?' T) + AVf T) c 

(15) 

where A is the Lagrange multiplier associated with the 
quality of the least-squares fit in the S zone. By aug- 
menting the cost function using a quantity that scales 
with the flux suppression of a potential companion, we 
introduce a damping term in the speckle suppression 
algorithm. This term forces the overall speckle reduction 
to compete with flux preservation in the S zone. Having 
chosen a positive definite cost function allows us to use 
standard quadratic programming routines to solve the 
following inverse problem: 

Find c DL and A DL that minimize J(c DL ,A DL ) 
under the constraint that c^f L ^ for all k 

where the subscript DL stands for "damped LOCI", 
d-LOCI hereafter. 

In practice we implement this optimization problem by 
nesting a linear least-squares solver within a one dimen- 
sional non-linear search, 



A 



DL 



miriA J(c(A), A) 



with: 



c(A) = (M^' fl) )t(v£f' T) 



A, 



(16) 



(17) 



where the superscript t stands for the linear least 
square minimization under the constraint that all coeffi- 
cients are positive. This algorithm will find the optimal 
weight between speckle and signal suppression and thus 
should behave better than regular LOCI with respect to 
the photometry and spectral features of a potential com- 
panion. Note that the damping in d-LOCI occurs in two 
separate instances: the positivity constraint on the coef- 
ficients and the conservation of the companion flux in the 
S zone. We used a set of synthetic companions to assess 
the relative performances of both LOCI and d-LOCI with 
respect to detectability. Our results are presented in Ap- 
pendix B, and illustrate that both methods yield compa- 
rable SNR. In the next section we delve into the problem 
of estimating of the spectro-photometry of companions 
detected in IFS data using LOCI type reductions. 

3. SPECTRAL EXTRACTION WITH D-LOCI 

3.1. Spectral bias and contamination 

If we write the discrete spectrum of a compan- 
ion, sampled at the spectral resolution of the IFS, as 
{Ip} p =i ... N\ i then the flux of a potential companion in 
a LOCI subtracted image at the wavelength X Po can be 
written as 



tC,LOCI _ tC 
Po Po 

jC,LOCI = jC 



L 



Bias 



Po 
Bias 



■ 1 Star — ICross (18) 
' 1 Star — c plpo,pIp • 

(19) 

pen 



The three biases in Eq. [19] are described as follows: 

• Ip Q zas corresponds to the amount of companion flux 
that is directly fitted by the PSF subtraction algo- 
rithm. It is the portion of the companion light 
that has been considered as starlight by the LOCI 
algorithm and thus subtracted. This bias can oc- 
cur in all implementations of LOCI, but is larger 
when using chromatic differential imaging, com- 
pared to ADI, because of the potential poor condi- 
tioning of the correlation matrix M^'^. The rela- 
tive size of the subtraction and optimization zones, 
which directly scales with Na, impacts Ip Q ias . This 
number ought to be carefully calibrated using syn- 
thetic point sources when measuring the photome- 
try within any LOCI processed image. 

• The bias Istar corresponds to the residual starlight 
after LOCI in the pixels of interest. If we assume 
that this residual background is constant around 
the companion PSF, it can be estimated and sub- 
tracted. As a consequence we neglect it here. 

• Icross corresponds to the contamination of the 
measured flux at a given wavelength by the com- 
panion flux present in the composite reference at 
adjacent wavelengths. Indeed the measurement at 
X Po is contaminated by the spectral content of the 
companion at nearby wavelengths Ip multiplied by 
the corresponding LOCI coefficient, c p , and a quan- 
tity, 7 Po ,p, that measures the spatial extend of the 
companion PSF between wavelength channels. 1Z is 
chosen such that for all p and q, j Pyq <C 1. However 
in the absence of damping the coefficients \c p \ can 



6 



Pueyo et.al 



become large enough to make the spectral contam- 



ination non negligible, | c P 7p 5< 



This spectral 



cross talk is specific to the use of LOCI with IFS 
images. It can be somewhat mitigated by optimiz- 
ing the value of N§ using synthetic companions. 
However such an optimization relies on knowing a 
priori the spectrum of the companion to charac- 
terize, and thus can only be carried out when this 
companion is bright enough to obtain a rough esti- 
mate of the spectrum in the absence of any speckle 
calibration. 



The final error on the retrieved spectrum is the sum of 
the algorithmic error, — I^ LOCI , and the spectro- 
scopic calibration uncertainties. Calibrating the P1640 
spectrograph involves a careful measurement of the in- 
strumental spectral response function, whose de t ermina - 
tion procedure is detailed in [Zimmer man et al. l (|2010a| ) ; 
Zimmerman et al.| ( |2010b| ); |Hinkley et al.| ( l2U10ap . The 
uncertainty on this response is not expected to be larger 
than 5 percents. As we will see below, the quantity 
Ip Q — I^ LOCI when using LOCI reduced data is much 
larger than this expected telluric calibration uncertainty 
and is thus very often the dominant source of errors. As 
a consequence we carried out our spectral fidelity compu- 
tations in "camera space", using wavelength series that 
have not been normalized by the instrumental spectral 
response function. By doing so we chose to leave aside 
discussions about telluric calibration and only focus on 
issues related to the effects of the PSF subtraction algo- 
rithms. 

3.2. PSF binary weighting 

In this section we consider the case of a faint compan- 
ion which has been detected and needs to be character- 
ized. We call V the image plane location of the detected 
off-axis source, and assume the geometry described on 
the left hand side of Fig. [I] and Fig. [2j where V and S 
have the same area but do not fully overlap. While the 
case of using as subtraction zone that is smaller than V 
is of particular interest for the detection problem, this 
configuration leads to large signal depletions and is thus 
not discussed in the present communication. The top 
panel of Fig. [I] illustrates a test case carried out using 
aperture photometry on an IFS data-set reduced usin g 
the LOCI algorithm described in |Lafreniere et a l. (2007). 
We introduced a fake companion located at 1" in 4 Al- 
cor data-cubes. This corresponds to about 10 minutes 
of exposure, each exposure being composed of N\ = 23 
frames (around 80 reference). We ran our reductions 
with constant LOCI parameters dr = 2, N$ = 1, g = 3 
and Na = 100. This choice of Na leads to solving in- 
verse problems which are slightly over-constrained and 
thus favorable to the regular LOCI in terms of detectabil- 
ity. In this figure the spectrum of the fake companion is 
similar to the one of the host star, except at 1.65/im, 
where a synthetic a 10% flux depletion was inserted, and 
its brightness was chosen such that its integrated delta 
magnitude over the J and H band was Amj + # = 6.5. 
The retrieved spectrum is very different from the injected 
one. Because the observed discrepancy cannot be de- 
scribed as a "gray gain" , it is a combination of hias and 
I cross i it cannot be calibrated using a set of synthetic 
companions. The 10% flux depletion at 1.65/im allows 



us to illustrate the impact of I cr0 ss- This weak absorp- 
tion is strongly amplified via wavelength cross-talk in the 
reduction method, as shown in the top panel of Fig. [I] 
This aspect is very difficult to calibrate using synthetic 
sources since this depletion depends on the depth of the 
putative absorption, which is unknown a-priori. We first 
sought to address these issues by exploring several binary 
weighting strategies of the LOCI cost function: 

• Subtraction zone masking consists on reducing I^as 
by not including the pixels in the S zone in the least 
square fit. It excludes the companion flux from the 
inverse problem and thus should preserve better the 
flux of potential companions. In practice it is im- 
plemented using a weighting function w(x,y) that 
is 1 in the optimization zone and in the subtrac- 
tion zone; and corresponds to solving the following 
optimization problem: 



Find c k st min 



o-s 



(20) 



The second panel of Fig. [T] illustrates the efficiency 
of this approach in our test case. I^as is indeed 
reduced, however it is achieved at the cost of an in- 
crease in I cr0 ss : the flux the 18 th spectral channel 
is strongly depleted because of the large coefficients 
at neighboring wavelengths. 

• Companion masking in the target image. One way 
to to reduce I cr0 ss consists on minimizing the mag- 
nitude of the Cfc's at the location of the compan- 
ion. Since the off-axis source has already been de- 
tected, this can be achieved by masking its pixels 
in the target image and thus solve the following 
least-squares problem: 



Find c k st min I \\w-p(x,y)T(x,y) 



k 



c k R k (x,y)\\ 



(21) 

where w-p(x, y) is equal to at the location of the 
detected companion and 1 everywhere else. This is 
equivalent to finding a composite reference which 
minimizes both the residual in the O zone, and the 
impact of the composite reference at the location 
of the companion: 



Find c k st min I \ \T — ^ c k R 



k\\0-V ' 



k\\V 



(22) 

This is illustrated on the third panel of Fig. [T] This 
approach can be interpreted as a simplified ver- 
sion of d-LOCI, without the coefficient positivity 
constraint and under the implicit assumption that 
A = 1. This approach considerably reduces Iuas 
to the point that the flux from the companion is 
overestimated. 

• Subtraction zone masking and companion masking 
in the target frame combines the two approaches 
and is illustrated on the second panel of Fig. [2] It 
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does alleviates spectral contamination but overes- 
timates the companion flux, albeit at a much lower 
level than solely masking the companion in the tar- 
get image. 

The analysis using synthetic companions presented on 
Fig. [I] and Fig. [5] was carried out at several angular 
separations and delta magnitudes and consistent results 
were obtained. While using subtraction zone masking 
in conjunction with companion masking in the target 
frame yields good spectral fidelity at low contrast lev- 
els, it degrades quickly when the synthetic companions 
gets fainter. Replacing the S zone masking by d-LOCI 
provides a robust alternative to this problem. 

3.3. Results 
3.3.1. Absolute and relative photometry 

IFS observations of faint companions with instruments 
such as PI 640 yield two astronomical observables: band 
averaged photometry and low resolution spectroscopy. 
Estimating both quantities from aggressively LOCI sub- 
tracted data-cubes require two fundamentally different 
reduction strategies: 

• Band-averaged photometry (J and H bands in the 
case of P1640): accurate broadband photometry 
is crucial for deriving effective temperatures and 
determining masses based on evolutionary models. 
The band averaged flux is insensitive to the intrin- 
sic spectral shape; this is regardless of whether 
or not the companion has absorption or emis- 
sion bands in the spectrum (main source of errors 
in Ibias)- Therefore the band-averaged LOCI re- 
sponse is relatively insensitive to underlying spec- 
trum and can be calibrated using synthetic com- 
panions with a flat spectra. This is equivalent t o 
work done by several authors | Marois et al. ( 2008b; 



iThalmann et^ j2009| > ; |Currie et ( |201 1| ) ; 



anson 



et al.| ( |201l[ ). LOCI parameters (7V a , dr\ his) can 
be tuned to optimize band- averaged performances. 
For instance, the bias of the integrated flux over 
both J and H bands, corresponding to the reduc- 
tions illustrated on Fig.[3J can be minimized using a 
LOCI parameters optimization based on synthetic 
companions. On this figure, regular LOCI reduc- 
tions (middle column of Fig. [31) exhibit more ac- 
curate band-averaged photometric estimates in the 
H band than in the J band: this suggests that the 
parameters can be further optimized in J band. 

• Low resolution spectroscopy: accurate spec- 
troscopy allows further constraints on physical pa- 
rameters such as effective temperature, surface 
gravity, and the detection of molecular absorption 
bands. Here the color response of the algorithm 
is sensitive to the intrinsic spectral shape (main 
source of errors in I cross)- Calibrating this re- 
sponse without foreknowledge of the companion's 
spectrum is difficult. This is the main problem ad- 
dressed in this paper: we show here that, even at 
the detection limit, d-LOCI provides an SED es- 
timates much more accurate than classical LOCI. 
Importantly, this does not require any calibration 
of the color response of the algorithm. 



We first discuss qualitatively the performance of four 
possibilities, with masking/no-masking and with regu- 
lar/damped LOCI, with respect to band-averaged pho- 
tometry and low resolution spectroscopy. To do so, we 
introduced a synthetic companion located at ~ 1" sep- 



aration in 4 d ata-cubes on the star Alcor, [Zimmerman 
et al. (2010b). This is about 10 minutes of 
with each 3D cube composed of N\ 



sure, 



total expo- 
= 23 wave- 
length frames. The 4 data-cubes allow a reference frame 
set of 80 frames. Unlike the A5V host Alcor, the 1.2 - 
1.8 micron SED of the synthetic companion has H-band 
excess; this lowers the companion signal-to- noise at short 
wavelengths (J-band) where the surrounding speckles are 
brighter. We then ran all reductions with constant LOCI 
parameters dr = 2, N§ = 1, g = 3 and Na = 100; 
this choice of Na leads to solving a set of inverse prob- 
lems that are slightly over-constrained, with 80 degrees 
of freedom and 100 constraints. The results are shown in 
Fig. [3J for three different companion magnitudes ranging 
a factor of 10 in flux. The leftmost column shows the 
before and after LOCI images used for these analyses. 
The "post-reduction" inset illustrates the companion's 
detection level in a single spectral channel. These images 
were obtained with classical LOCI and no priors on the 
companion's location. The behavior of the four different 
reductions illustrated in this figure allow insights into the 
algorithmic photometric and spectro-photometric accu- 
racy of both unmasked and masked approaches. We now 
refer the reader to Fig. [I] Fig. [2] and Fig. |3j 

Flux preservation properties oi unmasked 
LOCI: The extracted spectral profiles represented by the 
lines labelled "Companion non masked" in Fig. [3] corre- 
spond to the implementation of non-masked reductions 
with classical LOCI (center column), solution (a), top 
panel of Fig. [T] and with d-LOCI (right column), solution 
(d), top panelof Fig.[2j d-LOCI exhibits an improvement 
in the algorithmic throughput of a factor of two. As there 
is no prior about the companion's location, and as the 
SNR of both solutions are similar (§3), we find that d- 
LOCI preserves a larger portion of the companion's flux. 
Moreover the top set of panels of Fig. [^demonstrate that 
classical LOCI without binary weighting can lead to gray 
spectral transmitivities. This holds when t he companion 
is as bright as the speckles as is seen in Bo wler et al. 
(2010), where the authors used a set of synthetic com- 



panions to derive a gray algorithmic spectral response for 
HR8799 b. However this property quickly degrades with 
contrast and more robust methods are needed when the 
companion is fainter than the speckles. 

Spectro-photometric accuracy in masked or 
damped approaches: in the middle column of Fig. [3j 
the spectral profiles represented by the lines labelled 
"Companion masked" correspond to the implementation 
of LOCI with both subtraction zone masking and com- 
panion masking in the target frames. This refers to so- 
lution (e) in Fig. ^1 In the right panel the same spectra 
are generated by cFLOCI with companion maskingin the 
target frames. This refers to solution (f) in Fig. [2] The 
performance of both methods degrades as the compan- 
ion gets fainter. However, the former clearly results in 
larger spectro-photometric errors than the latter (with 
damping) . 

These findings are summarized in Table, [I] which 
describes the various levels of refinement of LOCI re- 
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ductions necessary to carry out accurate measurement 
of band-averaged photometry and low resolution spec- 
troscopy. We identified three levels of synthetic compan- 
ion brightness: brighter than speckles (can be detected 
in a broadband band image but its photometry is con- 
taminated by the speckles ), as bright as the speckles 
(can only be identified by its immobility in a movie of 
sequential wavelengths), fainter than the speckles. The 
level of refinement of LOCI reduction was qualitatively 
inferred from the series reductions carried out to gen- 
erate Fig. [3] (15 different companion brightness ranging 
from Am = 5 to Am = 8). Note that in this paper, 
no effort was made to optimize the sensitivity to faint 
off-axis sources. For details pertaining to the behavior 
of the detectability of PI 640 as a function of observing 
conditions (e.g seeing, host star magnitude, integration) 
and reduction procedure (e.g opti mization of LOCI pa - 



rameters), we direct the reader to Crepp et al. (2010) 



Next we quantify the absolute and relative spectral ac- 
curacies that are obtained using the two most promising 
approaches identified in Table. [I] dual PSF masking with 
classical LOCI and single PSF masking with d-LOCI. 

3.3.2. Band-averaged photometry absolute error 

Here the absolute spectro-photometric error is the root 
mean squared error in the measurement of the companion 
flux, integrated over the spectral channels spreading the 
J and H band. 



<^ Nx 1 (s telluric ) 2 (I c — Ip ,LOCI ) 2 
e Abs = — — — —7; (23) 



\^ N x Q telluricjC 
Z-^p=l *p 1 p 



where s 



telluric 



is the spectral response function, derived 



using ref erence A stars acc o rding to the method de- 
tailed in Zimmerman et al. (2010a). The left column 
of Fig. [4] shows how this metric varies as a function of 
the brightness and angular separation of a potential com- 
panion. The band-averaged photometric accuracy scales 
with the absolute spectro-photometric error. For com- 
parison, the average level of photometric uncertainty in 
current ADI broa dband discov e ries o f faint companions 



(5r 



0.2 



Marois et al. 



(2009); Currie efal 



LOCI consi 



(2008); 



Thalmann et al. 



anson et al.| (|201l|)) is in- 



dicated on Fig. [I] d-LOCl consistently leads to smaller 
absolute spectro-photometric errors for the reduction pa- 
rameters (Na, g, Ns,dr) tested here. The improvement 
in CAbs between dual masked LOCI and d-LOCI varies 
as a function of angular separation. This suggests, as 
discussed above, that the LOCI parameterization can be 
further optimized in order reduce the systematic error 
bar on band-averaged photometry. Note that most of 
the spectro-photometric error when using d-LOCI in the 
faintest case of Fig. [3] (bottom-right panel) consists of a 
band-average bias. Such an error can be mitigated using 
a set of synthetic companions to optimize LOCI param- 
eters. In the faintest case, the actual improvement in 
relative SED (up to a scaling constant) provided by d- 
LOCI, when compared to classical LOCI (bottom-middle 
panel of Fig. [3|, is actually significant. We quantify this 
improvement m the next sub-section. 

Note that the two reductions methods considered on 
Fig. [4] assume that a companion has already been de- 
tected and has been masked in the target frames. As a 



consequence they can only be carried out for sources al- 
ready identified using non-masked LOCI reductions. We 
separately evaluated the detectability of the synthetic 
companions in Fig. [4] and did not include sources below 
the band- averaged 5 — a level limit. We find that in the 
case of the P1640 data-set used for these calculations, 
this regime occurs at Am ^7.3 with slight variations 
depending on the angular separation. While this num- 
ber is relatively modest we emphasize that this is mostly 
due to our conservative choices of data-sets (average ob- 
serving conditions, no wavefront calibration before ob- 
servations) of fake spectrum (faint in the J band), short 
exposure times (9 minutes) and reduction parameters 
{Na,Q, Ns, dr). Further optimizations of theses parame- 
ters can surely lead to deeper detectability and absolute 
error that is consistent with the band averaged photo- 
metric error bars published in previous ADI discoveries; 
all the way to the 5 — a detection limit. 

3.3.3. Low resolution spectroscopy: relative error 

When seeking to characterize the spectral profile of a 
discovered companion within its near infrared bands, cal- 
ibrating the color response of the LOCI reduction with- 
out foreknowledge of the companion's spectrum is diffi- 
cult. When it is sufficiently bright with respect to the 
speckles, the special case of a gray spec tral response for 



Bowler et al. 



(2010) provides a thor- 



the reduction arises, 
ough illustration of this regime in the case of HR8799b 
When the algorithmic color response is gray, it is safe 
to assume that I^as is responsible for most of the flux 
depletion. Indeed, since I cr oss mostly amplifies local ab- 
sorption features, it is unlikely to be constant over the 
full spectrum. In this particular case both the photom- 
etry and the low resolution spectral features of the ob- 
ject can be estimated. However, as the companion gets 
fainter, the reduction's spectral response starts varying 
with wavelength. Without prior knowledge about the 
spectral profile of the companion, it is very difficult to 
disentangle whether such variations are the sole result of 
LOCI (e.g Ibias) or stem from the spectrum of the faint 
source propagated though LOCI (e.g Icross)- One solu- 
tion to address this issue would be to carry out a series of 
iterative reductions: (1) extract a biased spectrum from 
a preliminary LOCI reduction, (2) subtract this spec- 
trum from the non-reduced data, (3) run a secondary 
LOCI reduction, (4) iterate. Since the biases depend 
on the set of LOCI parameters chosen for the reduction 
(Na, g, N$,dr), this process should be repeated over a 
wide set of parameters in order to ensure that the biases 
on the SED are minimized. In this paper we answer an 
orthogonal question. Namely, how to minimize the rela- 
tive spectro-photometric bias in one iteration, for a well 
chosen but not optimized set of LOCI parameters? 

We quantify the relative spectro-photometric error as 
measure of the uncertainty in estimating the Spectral 
Energy Distribution of the discovered object: 



tRel 



Nx^2(s t p elluric ) 2 



T C,LOCI 



IS 



jC \r^N x Q telluricjC 
1 k l^k=l b p 1 ' 

(24) 



Q telluric 

k=l *p x k l^k=\ °p x k 



In other words €r 6 i captures of how well the overall shape 
of the low resolution spectrum can be retrieved, up to a 
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scaling constant. Assuming that the band- averaged pho- 
tometry of the object has already been estimated, this 
SED provides the quantitative basis to further charac- 
terize the discovered companion. In particular it pro- 
vides remarkable level arm to estimate both surface grav- 
ity, atmospheric c hemis try |Burgasser et al.| ( |2QQ6[ ) ; |Op-| 
penheimer et al. (1995) and, to a lesser exte nt, the at- 
mospheric c loud content of the faint object |Burgasse"r| 
et al. (2010). Therefore the spectro-photometric error 
will have a direct impact on the estimation of the physical 
parameters of the faint exo-planets discovered by cam- 
paigns such as GPI, SPHERE, P3000-P1640. Without 
any priors about the spectral content of such objects, cal- 
ibration of the biases introduced by the reduction using 
synthetic sources becomes a very delicate exercise. A first 
avenue that can be explored when seeking to address this 
issue of faint companion spectral characterization is to 
devise a new class of al gorithms that rely on prior s about 



campaigns JGraham etH] (|2007|); iDo hle n et al] ( |2006D ; 
Dekany et al.| ( |200(j[ ); |Hinkley et ai.| ( |20rT| ~ 



the speckle field (e.g. Burke & Devaney (2010); 



Mesa 



et al. ( 2011[ )) or the spectra l type of the companio n (e.g. 



spectral deconvolution as in Sparks & Ford ( 2002[) ). The 
solution proposed in the present manuscript, d-LOCI, fo- 
cuses class of solution which assumes no priors about the 
quasi-static errors or the companion. 

The right panel of Fig. [4] demonstrates that d-LOCI 
is a robust approach when seeking to characterize the 
spectral features of a faint companion without any fore- 
knowledge of its spectral type. Indeed, the relative error 
using d-LOCI in conjunction with target frame mask- 
ing is, for most separations and levels of contrast, half 
an order of magnitude smaller than the one obtained 
using dual masked classical LOCI. More importantly it 
yields SED uncertainties below the error-bar level of state 
of the art ban d-averaged photom etric es timates (|Marois1 
et al.| ( |2008D ; |Thalmann et al.| fl2009|) ; |Metchev eOT] 



2009 ) ) . This feature is quite remarkable since it re- 
uces the spectro-photometric uncertainties associated 
with observations with an IFS without field rotations to 
the level of observations carried out in a more favorable 
mode: ADI (I cr oss = 0) and images integrated over sev- 
eral P1640 channels (more photons from the companion). 
For the brightest synthetic companions tested here, the 
SED bias obtained using d-LOCI is of the order of the 
P1640 telluric calibration uncertainties. We also empha- 
size that such results are obtained without any iterations 
on LOCI parameter or underlying companion spectrum, 
and thus d-LOCI reduced images are a good first guess 
when seeking to obtain even higher spectro-photometric 
accuracy using iterative schemes. 

d-LOCI is thus an efficient method when one seeks to 
estimate the SED of a detected point source in P1640 
data, without any prior knowledge about the spectral 
profile of the companion or about the chromaticity of 
the post coronagraphic speckles. This result ought to be 
carefully extrapolated to other instruments, in particular 
using similar tests based on synthetic companions. We 
expect however that such conclusions will stand when 
the raw contrast of direct imaging campaigns wil l im 



prove, either usi ng a real-time calibration system i Wal- 
lace et al. (2009)) or advanced pre-observing calibration 
methods (Serabyn et al. ( 2010[ )). We thus envision that 
the approach suggested in herein will be critical to the 
estimation of the physical parameters and atmospheric 
content of exo-planets detected via future direct imaging 



4. CONCLUSION 

In this paper we discuss an algorithm for Point Spread 
Function subtraction in data acquired with an Integral 
Field Spectrograph, that is efficient and accurate for 
retrieving the spectral content of faint companions lo- 
cated within several diffractio n widths of a target star. 
A previous communication by I Crepp et al. (|20 1 [) used 
the newly fu nctional Project 1640 spectrograph,"^* 



|et al.| ( |2011[ ) , as a laboratory to devise a robust met hoc 
purely tor off-axis faint source detection. This method 
relies on incorporating the spectral diversity inherent in 
IFS data into a Locally Optimized Combination of Im- 
ages (LOCI) PSF calibration scheme. In that paper, we 
reported H-band 5 a detection limits of ~ 10 -5 in con- 
trast, at 1 arcsec separation, in a 20 min of exposure on 
an H = 3 star. In the present paper, we expand upon the 
general LOCI scheme, and discuss the spectral character- 
ization of already detected faint companions. We show 
that when LOCI is tuned to generate an aggressive com- 
posite calibration PSF, the spectro-photometric content 
of the faint object is biased. It is affected by both signal 
depletion, which is common to all LOCI implementa- 
tions, and by spectral cross-talk with neighboring wave- 
lengths, which is specific to IFS imaging spectroscopy. 
Herein, we introduced a new algorithm, "damped LOCI" 
(d-LOCI) that alleviates many of these shortcomings. 

We modified the formalism underlying LOCI PSF cal- 
ibration by introducing an augmented cost function that 
seeks to preserve the flux from the companion in ques- 
tion. We compared the efficiency of LOCI and d-LOCI 
for "detectability" , and found that both approaches per- 
formed similarly in the various test cases. We identified 
second order trends in the behavior of SNR as a function 
of exposure time and LOCI parameters. Such trends 
can be explained by the conditioning of the inverse prob- 
lem that is at the heart of LOCI reduction, and have 
only marginal bearing on t he conclusions perta ining to 
"detectability" discussed in Crepp et al.| p010[ ). While 
d-LOCI clearly increases the signal ol taint companions, 
it results in a noise floor that is shallower than that ob- 
tained in LOCI, and thus the detectability metric of both 
approaches is similar. 

We finally delved into the problem of spectral extrac- 
tion. We first used synthetic companions to explore the 
behavior of classical LOCI reductions, implemented with 
a variety of binary weightings. Our tests showed that a 
classical imple mentation of LOCI, w ith well chosen pa- 
rameters as in |Lafreniere et al" (2007), can be applied to 
companions slightly brighter than the speckles and yields 
an "almost gray" algorithmic throughput that yields a 10 
percent error in SED estimates. In this regime we find 
the spectral fidelity can be improved by the use of a dual 
binary weighting when posing the PSF reconstruction 
problem. Since we expect telluric calibration errors be- 
low 10 percent, the main source of errors in the spectro- 
photometric characterization of the discovered compan- 
ions resides in bias within the reduction algorithm. More- 
over these biases increase rapidly as function of contrast, 
and when the companion is fainter than the speckles, the 
error in the SED estimate can be as large as 50 percents. 
Using d-LOCI, in conjunction with a binary weighting 
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of the target frames at the location of the companion, 
consistently reduces this error by a factor of five: it low- 
ers the algorithmic errors in the SED close to the levels 
imposed by the telluric calibration. Translating these 
uncertainties in error bars on the estimates of the phys- 
ical parameters of potentially discovered companions is 
beyond the scope of this paper, and will be illustrated in 
a future communication. All results are obtained using 
the Project 1640 IFS data, but we fully expect that these 
results can be generalized to future high contrast imag- 
ing instruments operating in conjunction with advanced 
wave front calibra tion techniques. These next generation 
instruments (IGraham et al.| (|2007|); | Dohlen et aL|j2006[ ); 

wnFol 



nqi 
J(g_ 

Dekany et at] ( |2UU6| ); jHinkley et al.| ( |2Ull| )) will oper- 



ate at contrasts significantly higher than those possible 
with the current PAL AO+P 1640 configuration, and they 
will use Integral Field Spectrographs as science cameras. 
Accurate in-situ spectral characterization of faint exo- 
planets buried in the residual speckle fields of the long 
exposures will thus be a crucial component of upcoming 
campaigns. It will lead to far tighter constraints on plan- 
etary bulk physical properties and atmospheric chemistry 
and would be a direct application of the d-LOCI methods 
discussed herein. 
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(a) Regular LOCI: 

Minimize least square fit in an optimization zone that includes the subtraction zone. 
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1.1 1.2 1.3 1.4 1.5 1.6 1.7 

(b) S-zone masking LOCI: wavelength 

Minimize least square fit in an optimization zone that does not includes the subtraction zone. 




1.1 1.2 1.3 1.4 1.5 1.6 1.7 
Wavelength 

(c) Companion masking LOCI: 

Minimize least square fit of the target frame in an optimization zone that includes the 
subtraction zone but not the pixels of the companion. 




?.1 1.2 1.3 1.4 1.5 1.6 1.7 
Legend: Wavelength 



— ^ Spectrum of the star, similar to the spectrum of the injected synthetic companion 

% " except for synthetic absorbtion at 1 .65 microns 
™^ Spectrum of the injected companion 

~ e ~ Sla Extracted Spectrum after LOCI 

Fig. 1. — Illustration of three PSF binary weighting approaches, 
with their respective performances in the case of a companion 
slightly fainter than the P1640 speckles, Am = 6. Note that the 
scale of the bottom panel has been compresses to accommodate 
for the large bias. None of these solution provide robust spectro- 
photometric accuracy: they all lead to combination of Ibias an d 
Icross- O is the optimizations zone, used to find the best least- 
squares fit of the speckles realization. S is the subtraction zone. V 
is the position of the planet. 7Z is the ensemble of reference images. 
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(d) d-LOCI 

Minimize least square fit of the target frame in an optimization zone while maximizing the 
residual in the subtraction zone. 
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(e) Companion and S-zone masking LOCI: 

Minimize least square fit of the target frame in an optimization zone that excludes the 
subtraction zone and the pixels of the companion. 
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(f) Companion masking d-LOCI: 

Minimize least square fit of the target frame in an optimization zone which excludes the pixels 
of the companion, while maximizing the residual in the subtraction zone. 
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Legend: 1.1 1.2 1.3 1.4 1.5 

Wavelength 

— ^ Spectrum of the star, similar to the spectrum of the injected synthetic companion 
% * except for synthetic absorbtion at 1 .65 microns 
Spectrum of the injected companion 

" e ~ e ** Extracted Spectrum after LOCI 

Fig. 2. — Illustration of three PSF binary weighting approaches 
combined with d-LOCI, with their respective performances in the 
case of a companion slightly fainter than the P1640 speckles, Am = 
6. Solutions (e) and (f) provide the best two spectral estimates and 
their biases are carefully quantified in the present manuscript. O is 
the optimizations zone, used to find the best least-squares fit of the 
speckles realization. S is the subtraction zone. V is the position 
of the planet. 1Z is the ensemble of reference images. The upwards 
arrows illustrate the flux preservation constraint in d-LOCI. 
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Fig. 3. — Qualitative comparison of the spectro-photometric ac- 
curacy of different LOCI reductions. We used synthetic compan- 
ions at 1 arcsec at three different levels of contrast. The injected 
spectrum is represented by a solid line in all panels. Left Panel: 
Images before and after LOCI. Middle panel: Classical LOCI, 
solution (a) in Fig. ^ and Fig. p] in crossed-dashed line; classi- 
cal LOCI with dual masking, solution (e), in circled dashed line. 
Right panel: d-LOCI, solution (d), crossed-dashed line, d-LOCI 
with target frame masking, solution (f), in circled dashed line. 
When using d-LOCI, solution (f), most of the residual error cor- 
responds to a gray gain, or absolute error. The relative SED error 
is smaller in this d-LOCI than in the classical masked case, solu- 
tion (d). This figure illustrates approaches which do not include 
any iterations on the underlying companion's spectrum or LOCI 
parameters 



14 



CO 
O 



-0.5 



log(Absolute Error) 

iog(B Abs ) = f(Am) 

Companion at 0.4" 




CD 
O 



-0.5 



Pueyo et.al 

log(Relative Error) 

log(E Re| )=f(Am) 
Companion at 0.4" 




0.4 Companion at 1 ' 

0.2 



CO 

o 




CO 

-0.5 



5 5.5 

Companion 




-0.5 



Companion at 1.5" 




A m 



A m 



O o Damped LOCI with companion masking in the target image 

x x Regular LOCI with companion masking in the target image and S zone masking 

Fig. 4. — Quantitative comparison of the spectro-photometric 
accuracy for LOCI reductions (d) and (e) as a function of contrast 
at three different angular separations. Left: band-averaged photo- 
metric error. Right: bias in the SED estimate. The gray area cor- 
responds to published band-averaged photometric error bars on re- 
cently discovered faint companions using ADI: 5m j &5mH < ±0.2 
mags. Ev en in the challenging "IFS +no field rotation" configura- 
tion, d-LOCI delivers systematic errors on the SED of the order 
of state of the art uncertainties derived from band- averaged ADI 
observations, this all the way to the detection limit of the PI64O 
instrument. 
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TABLE 1 
SUMMARY OF OUR FINDINGS 



Preferred LOCI method ( a )for given companion brightness and science goals 


Companion brightness 


Photometry w/o iterations^ Photometry with iterations^ SED estimates'^ 


Brighter than speckles 


dual masked LOCI LOCI dual masked LOCI 


As bright as speckles 


d-LOCI dual masked LOCI d-LOCI 


Fainter than speckles 


NONE d-LOCI d-LOCI 



Note: (a)LOCI reduction with the smallest degree of sophistica- 
tion that delivers a given science goal for a given companion 
brightness. 

(b) The iterations can be carried out on the LOCI parame- 
ters or the brightness of the object. In any case a sequence 
of synthetic companions injections and LOCI reductions is 
necessary in this scheme. 

(c) e.g one single LOCI reduction is sufficient to provide a 
robust photometric estimate. 

(d) Here we assume that LOCI biases on the SED depend 
on the underlying spectrum, and thus are very difficult to 
calibrate using iterations without any prior on the faint 
companion. 
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APPENDIX 



APPENDIX: FRAME SELECTION AND GEOMETRY 

In this communication we consider the case of a series of IFS cubes that have already been centered and spatially 
re-scaled so that the optical a r tifact s are stationary within the image space and the astronomical image has been 
stretched radially, Crepp et al. (2010); Sparks & Ford (2002). Proper scaling and centering of the images is crucial to 
subsequent data analysis. We carry out this stage by finding the centering and scaling parameters that maximize the 
auto-correlat ion between two images. We solve this non-linear search using the sub-pixel image registration algorithm 
described by |Guizar- Sicairos et al. l (120 08), which has been modified using the Discrete Fourier Transform algorithm 
described in Soummer et al. ( |200yd| ), in order to accommodate for the scaling parameter. The series of images 
resulting from such a calibration corresponds to a virtual four dimensional cube, with two spatial dimensions, one 
spectral dimension and a temporal dimension. We seek to devise a reduction algorithm that will stack the signal of 
a potential companion for a given wavelength, while subtracting the artifacts due to optical imperfections. We write 

where the subscript p = 1 ... N\ stands for the spectral channel, and 

Associated with these slices 

is a collection of scaling parameters addressed similarly, jp Q \ We first build an integrated target frame T, with an 
associated scaling parameter j Po , that corresponds to the sum of the frames at the wavelength of interest A, 



each frame in our four dimensional cube as Ip Q 
the superscript q = l...N cu b es stands for the cube number, or the temporal dimension. 



Po- 



N cubes 

T= E 4? 

q=l 

1 N cubei 

iy cubes 



<v(9) 
'Po 



(Al) 
(A2) 



This is illustrated in the top panel of Fig. |5j where slices of the same wavelength are co-added. We seek to create 
a composite reference that locally subtracts the quasi-static speckles without removing the flux from a putative 
companion. Thus, at a given radial location within the target frame, written as r, we need to sort the remaining 
frames into two subsets: {Rk}, the subset of slices that will be used to create the composite reference and {J/e}, a 
subset of "useless" frames whose wavelength is too close to X p so that the image of a potential companion would 
overlap with the region of interest. Since the several wavelengths in the data-cube provide radial diversity, the image 
of a potential companion m oves radially in the rescaled space, we can rewrite the self-subtraction condition established 
by Lafreniere et al. (2007) as: 



K = {R k } = {#> : |r( 7 <*> - J Po )\ > N S W} 
J = Uk} = {I^ ■■ |r( 7 ^ - 7p„)I < N S W} 



(A3) 
(A4) 



where W is the Full Width Half Maximum (FWHM) of a PSF at wavelength A Po , and N§ is a parameter that determines 
the aggressiveness of the LOCI reduction algorithm. This sorting process involves rejecting the slices whose wavelength 
is too close to X Po from the collection of references, so that there is minimal companion flux in the composite references. 
We actually do not integrate over time the slices in the reference ensemble and we let LOCI choose the best coefficients 
to form an optimal reference PSF over time and wavelength. The strategy summarized on Fig. [5] was chosen on order 
to mitigate computational cost. Indeed, calculating a LOCI reference PSF for each slice of each exposure using PSFs 
at other wavelengths and slightly de-correlated in time, would mean running N cu ij es x N\ LOCI reductions, with each 
reduction looking for coefficients over an ensemble of N cu b es x (1 — e(N§))N\ reference PSFs (c(N§) is the number 
of slices excluded to avoid cross talk). Instead our reduction strategy only leads to N\ LOCI reductions with each 
reduction looking for coefficients over an ensemble of N cu b es x (1 — e(Ns))N\ reference PSFs. The performances of 
each reductions strategy in terms of detectability depends on the actual shot noise level of the speckle field (e.g. 
whether LOCI actually reconstructs the true realization of the speckle field in the target image, or fits the photon 
noise associated with it). The true merit of a given frame selection strategy thus depends strongly on observational 
conditions such as stellar magnitude or integration time. Finding the optimal binning scheme for given observational 
conditions is beyond the scope of the present manuscript and we did not delve into such trade-offs. Instead, we chose 
a hybrid frame selection approach which preserves the temporal information in the ensemble of reference PSFs while 
integrating the time variability of the target science image. The next step of the reduction algorithm is to find a series 
of optimal coefficients in order to build the composite reference. This is the stage that can lead to a loss of the spectral 
information of the detected companion. In order to formally de scribe the challenges of LOCI reductions we first briefly 
review the notations involved in the framework introduced by Lafreniere et al. (2007). As shown on Fig. 2, the local 
composite reference PSF is constructed by finding a linear combination ol reference frames, {Rk G 7^}, that yields 
the best least-squares fit of the speckle pattern in the target image, T. This locally optimal calibration image is then 
subtracted from the target image. However, in order for this reference to only describe the speckle field and not the 
flux of a potential companion, its coefficients ought to be determined in a zone of the image that is larger than the 
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zone within which the subtraction will occur. Thus the LOCI algorithm uses two regions of interest in the images, a 
subtraction zone, <S, where we will subtract the composite reference from the target, and a large optimization zone, 
0, that surrounds <S, t hat we will use to find the weights of the composite reference. We follow the presentation of 
Lafreniere et al. (2007) and use an annular geometry, as illustrated in Fig. [6] The optimization zone is defined by its 
area, expressed m number of PSF cores, Na, and its aspect ratio g. The aspect ratio is a unit-less quantity that scales 
with the radial elongation of the zone (9, 

Ar 

9 = (r + Ar/2)A0 (A5) 

where r is the radial location of the inner annulus of the optimization zone, Ar its radial extent and A</> its azimuthal 
extent. They can be written as a function of Na, 9 and W: 



Ar = 
Acj)-- 



7rN A gW 2 



4 
2r 

W y irN A 



1/2 



(A6) 
(A7) 



In the geometry chosen here, illustrated in Fig. [6| the inner radius of the subtraction and optimization zones is similar 
and equal to r. Their azimuthal span is likewise identical and equal to A(f>. The radial extend of S zone is dr while 
the O zone stretches all the wa y to r + Ar. The efficien cy of the LOCI reduction strongly depends on the choice of the 
parameters N A , dr and N$. Lafreniere et al. (2007) explored the parameter space for LOCI applied to ADI, and 
Crepp et al. (2010) presented similar work in the case of chromatic differential imaging. In this paper we show that 
a standard LOCI implementation applied to IFS data yields very good results in term of detectability, but perturbs 
the spectrum of the companions discovered. While the errors associated with the spectral estimation can be reduced 
by exploring the N A ,g,dr,Ns parameter space, we propose a more direct solution that consists of constraining the 
problem associated with the LOCI coefficient search. 



APPENDIX: DETECTION AND SIGNAL TO NOISE RATIO 

We first study the impact of d-LOCI on the detectability of faint point sources and show that the damping 
introduced herein only marginally affects the SNR of detected companions. To do so, we use a data set obtained on 
the bright A5V star (V = 4) Alcor on March 16 2009. We introduced in the data a spiral of 11 synthetic companions 
with a flat spectrum, at angular separations varying from 0.3 to 0.5 arcse c. The inject i on of the artificial sources 
occurs as high up as possible in the cube processing pipeline, as described in Crepp et al. ( 20101). Alcor br i ght vi sible 
magnitude enables good performances from the Palomar Adaptive Optics system (PALAO, [Dekany et al. 1 (2006)). A 
faint M-dwarf was recently detected using the P1640 spectrograph, |Zimmerman et a~ ( 2010bp , maki ng th is data set 



an ide al test case. The contrast of the astronomical companion is moderate, 1/30 as estimated by |Mamajek et al. 



(2010), and the secondary star is visible above the speckle halo: it provides a second reference for the photometry 
ol the synthetic companions inserted. Fig. 7] presents unprocessed images and the results of a data reduction, both 
integrated in J and H-bands. Since the speckles are brighter at shorter wavelengths, and the spectrum of the synthetic 
sources injected is gray, the overall contrast associated with them varies between 10 -2 and 10 -3 in the H band, and 
10 -3 and 10 -4 in the J band. As a consequence some companions are too faint to be detected in the J reduced images. 
These results were obtained using the the following parameters: one exposure of 140 seconds, N\ = 23 slices across 
the cube of interest, N A = 100, g = 1, dr — 2 and N$ = 1. This choice of parameters, in particular g = 1, implies that 
two adjacent synthetic companions might be located in the same LOCI zones, and thus might artificially hamper the 
SNR of reduced images. However we included our "conservative", g = 1, results in this manuscript, as its purpose 
was not do quantify the ultimate detection limits of PI 640. We refer to Crepp et al. ( 2010| ) for an example of a PI 640 
reduction when introducing the companions one at a time. Note also that the exposure time can be lengthened and 
the other algorithm parameters can be finely tuned for better SNR. However we do not delve into th e exploration of 
the whole parameter space: the results of such an optimization are presented in Crepp et al. (2010). In this section 
we describe a couple of test cases which illustrate the different regimes of the regular and d-LOCl algorithms with 
regard to point source detectability. 



We first consider the case of an over-constrained problem for which the number of free parameters, which scales 
with the number of reference frames, is smaller than the number of constraints, which scales with Na, the number 
of PSF cores in the zone O. We only reduce one data cube. This corresponds to a short integration, 140 seconds, 
and a small number of references to choose from in order to build the calibration PSF: N-ji < N\ = 23 (the actual 
value of N-jz varies as a function of angular separation, at the location of the injected synthetic sources a typical 
value is N-r ~ 15). We fix the geometry of the problem as follows: g = 1, dr = 2 and 5 = 1. We then vary the size 
of the optimization zone by changing the value of N A from 20 to 300, and thus study the behavior of both LOCI 
and d-LOCI in the over-constrained case since N A > Nn in these configurations. Fig. [8] illustrates how the average 
SNR of the 11 injected companions, integrated over the wavelength channels in both J and H bands, behaves in this 
configuration. The signal is estimated using aperture photometry at the known location of each artificial companion. 
The noise is estimated by computing the standard variation of the flux of a set of 100 boxes surrounding the spiral 
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of fake planets. Each photometric box has the area of a PSF core. Since the purpose of this section is to illustrate 
general trends pertaining to detectability with LOCI and d-LOCI, we chose to only isolate two "raw contrast" levels: 
"faint" companions in J band and "bright" companions in H band, and only show how the averaged SNR over 
the 11 companions behaves of for these two sub-groups. A detailed study, with synthetic sources of continuously 



decreasing brightness, was carried out in |Crepp et al. ( |2Q1Q[ ) in order to derive the sensitivity limit of P1640. Both 
LOCI and d-LOCI present the same overall level ol performance with slight upward SNR trend in both bands. In 
spite of having a penalty term that worsens the least-squares residual in the zone of the image where the PSF is 
subtracted, d-LOCI performs nearly as well than regular LOCI in the over-constrained case. The small SNR difference 
can be understood as follows: the noise suppression using pseudo-inverse methods is always more aggressive than 
when a damping term in introduced, and, in the over-constrained case, the likelihood of fitting companion flux with 
starlight is low. The upward trend in SNR as a function of the area of the optimization zone means, for the data-set 
studied here, that the higher companion throughput (e.g. larger signal) achieved by using a larger Na slightly dom- 
inates the degradation of the noise due the fact that the composite calibration image is calculated over a less local area. 

We then study the case of an under-constrained problem, for which the number of reference frames is larger than 
the number of PSF cores in the zone O. To do so, we now fix the area of the optimization zone Na = 100 and vary 
the number of cubes reduced from 1 to 14. According to our frame selection strategy, illustrated on Fig. |5j using 
more data-cubes should increase the signal, since co-adding cubes increases the exposure time, and should also reduce 
the noise after LOCI, since Nn increases with the number of cubes, thus making the subtraction more efficient. This 
means that we expect the SNR to be an increasing function of N-ji- Our result are shown on Fig. [9j Classical LOCI 
H-band results exhibit this behavior, with SNR increasing from 8 to 11, and yield better performances than d-LOCI 
reductions, for which the SNR remains constant at 8. Note that the rise in detectability does not scale as y/Nq^. 
This suggest that both our frame selection method (that was chosen for computational cost purposes) our LOCI 
parameters are not fully optimized for the synthetic companions injected. In the J-Band, for which the contrast 
is less favorable, d-LOCI leads to constant detectability levels (SNR = 4.5), while the performances of classical 
LOCI actually degrade with Njz- In this configuration the likelihood of fitting companion flux is high, since there 
are more degrees of freedom than there are constraints, and thus the algorithmic throughput decreases when N<jz 
increases. When the companion is faint enough (in the J-band) then classical LOCI implementations that are under- 
constrained can reach a regime for which the companion's flux suppression is greater than the speckle noise attenuation. 

The cases illustrated on Fig. [5] and Fig. [9] first confirm that when seeking for faint companions using LOCI reduction 
methods one should optimize the algorithm parameters beforehand using synthetic sources. In particular they show 

(R R) 

how the conditioning of the correlation matrix, M^, ' ; is a critical metric in this parameter search and that Na and 
Nft ought to be adjust e d so that the inverse problem is not under-constrained. This was carried out in the case of 



P1640 by Crepp et al.| (2010). Second, they show that d-LOCI detectability performances are close to the ones of 
classical LOCI. Overall, the damping makes d-LOCI less sensitive to the conditioning of the inverse problem, and his 
flat response to variations of the algorithm parameters might be used in order to discriminate astronomical companions 
from speckles residuals in the case of marginal detection. However a careful study of this second order trend is beyond 
the scope of the present paper. 

ACKNOWLEDGEMENTS 

We thank our anonymous referee for his/her very useful suggestions and Christian Marois for very productive 
conversations during the early stages of this manuscript. The research described in this publication was carried out 
in part at the Jet Propulsion Laboratory, California Institute of Technology, under a contract with the National 
Aeronautics and Space Administration. Project 1640 is funded by National Science Foundation grants AST-0520822, 
AST-0804417, and AST-0908484. This work was partially funded through the NASA ROSES Origins of Solar Systems 
Grant NMO710830/102190, the NSF AST-0908497 Grant. The adaptive optics program at Palomar is supported 
by NSF grants AST-0619922 and AST-1007046. LP was supported by an appointment to the NASA Postdoctoral 
Program at the JPL, Caltech, administered by Oak Ridge Associated Universities through a contract with NASA. LP 
and SH acknowledge support from the Carl Sagan Fellowship Program. 



REFERENCES 



Aime, C, & Soummer, R. 2004, ApJ, 612, L85 

Biller, B. A., Close, L. M., Masciadri, E., Nielsen, E., Lenzen, R., 

Brandner, W., McCarthy, D., Hartung, M., Kellner, S., 

Mamajek, E., Henning, T., Miller, D., Kenworthy, M., & 

Kulesa, C. 2007, ApJS, 173, 143 
Borde, P. J., & Traub, W. A. 2006, ApJ, 638, 488 
Bowler, B. P., Liu, M. C, Dupuy, T. J., & Cushing, M. C. 2010, 

ArXiv e-prints 

Burgasser, A. J., Burrows, A., Sz Kirkpatrick, J. D. 2006, ApJ, 
639, 1095 

Burgasser, A. J., Simcoe, R. A., Bochanski, J. J., Saumon, D., 
Mamajek, E. E., Cushing, M. C, Marley, M. S., McMurtry, C, 
Pipher, J. L., & Forrest, W. J. 2010, ApJ, 725, 1405 

Burke, D., & Devaney, N. 2010, J. Opt. Soc. Am. A, 27, A246 



Crepp, J. R. and Pueyo, L. and Brenner, D. and Oppenheimer, 
B. R. and Zimmerman, N. and Hinkley, S. and Parry, I. and 
King, D. and Vasisht, G. and Beichman, C. and Hillenbrand, L. 
and Dekany, R. and Shao, M. and Burruss, R. and Roberts, 
L. C. and Bouchez, A. and Roberts, J. and Soummer, R.. 2011, 
ApJ, 729, 132-+ 

Currie, T., Burrows, A., Itoh, Y., Matsumura, S., Fukagawa, M., 
Apai, D., Madhusudhan, N., Hinz, P. M., Rodigas, T. J., 
Kasper, M., Pyo, T., & Ogino, S. 2011, ApJ, 729, 128 

Dekany, R., Bouchez, A., Britton, M., Velur, V., Troy, M., 
Shelton, J. C, & Roberts, J. 2006 in (SPIE), 62720G 



IFS reduction, Spectral Extraction 



19 



Dohlen, K., Beuzit, J.-L., Feldt, M., Mouillet, D., Puget, P., 
Antichi, J., Baruffolo, A., Baudoz, P., Berton, A., Boccaletti, 
A., Carbillet, M., Charton, J., Claudi, R., Downing, M., 
Fabron, C, Feautrier, P., Fedrigo, E., Fusco, T., Gach, J.-L., 
Gratton, R., Hubin, N., Kasper, M., Langlois, M., Longmore, 
A., Moutou, C, Petit, C, Pragt, J., Rabou, P., Rousset, G., 
Saisse, M., Schmid, H.-M., Stadler, E., Stamm, D., Turatto, M., 
Waters, R., & Wildi, F. 2006, in Presented at the Society of 
Photo-Optical Instrumentation Engineers (SPIE) Conference, 
Vol. 6269, Ground-based and Airborne Instrumentation for 
Astronomy. Edited by McLean, Ian S.; Iye, Masanori. 
Proceedings of the SPIE, Volume 6269, pp. 62690Q (2006). 

Fitzgerald, M. P., & Graham, J. R. 2006, ApJ, 637, 541 

Give'on, A., Belikov, R., Shaklan, S., & Kasdin, J. 2007, Opt. 
Express, 15, 12338 

Graham, J. R., Macintosh, B., Doyon, R., Gavel, D., Larkin, J., 
Levine, M., Oppenheimer, B., Palmer, D., Saddlemyer, L., 
Sivaramakrishnan, A., Veran, J. -P., & Wallace, K. 2007, ArXiv 
e-prints, 704 

Guizar-Sicairos, M., Thurman, S. T., & Fienup, J. R. 2008, Opt. 
Lett., 33, 156 

Guyon, O., Pluzhnik, E., Martinache, F., Totems, J., Tanaka, S., 
Matsuo, T., Blain, C., & Belikov, R. 2010, PASP, 122, 71 

Hinkley, S., Oppenheimer, B. R., Brenner, D., Zimmerman, N., 
Roberts, L. C., Parry, I. R., Soummer, R., Sivaramakrishnan, 
A., Simon, M., Perrin, M. D., King, D. L., Lloyd, J. P., 
Bouchez, A., Roberts, J. E., Dekany, R., Beichman, C., 
Hillenbrand, L., Burruss, R., Shao, M., & Vasisht, G. 2010a, 
ApJ, 712, 421 

Hinkley, S., Oppenheimer, B. R., Soummer, R., 

Sivaramakrishnan, A., Roberts, Jr., L. C., Kuhn, J., Makidon, 
R. B., Perrin, M. D., Lloyd, J. P., Kratter, K., & Brenner, D. 
2007, ApJ, 654, 633 

Hinkley, S., Oppenheimer, B. R., Zimmerman, N., Brenner, D., 
Parry, I. R., Crepp, J. R., Vasisht, G., Ligon, E., King, D., 
Soummer, R., Sivaramakrishnan, A., Beichman, C., Shao, M., 
Roberts, L. C., Bouchez, A., Dekany, R., Pueyo, L., Roberts, 
J. E., Lockhart, T., Zhai, C., Shelton, C., & Burruss, R. 2011, 
PASP, 123, 74 

Janson, M., Carson, J., Thalmann, C, McElwain, M. W., Goto, 
M., Crepp, J., Wisniewski, J., Abe, L., Brandner, W., Burrows, 
A., Egner, S., Feldt, M., Grady, C. A., Golota, T., Guyon, O., 
Hashimoto, J., Hayano, Y., Hayashi, M., Hayashi, S., Henning, 
T., Hodapp, K. W., Ishii, M., Iye, M., Kandori, R., Knapp, 
G. R., Kudo, T., Kusakabe, N., Kuzuhara, M., Matsuo, T., 
Mayama, S., Miyama, S., Morino, J., Moro-Martm, A., 
Nishimura, T., Pyo, T., Serabyn, E., Suto, H., Suzuki, R., 
Takami, M., Takato, N., Terada, H., Tofflemire, B., Tomono, 
D., Turner, E. L., Watanabe, M., Yamada, T., Takami, H., 
Usuda, T., & Tamura, M. 2011, ApJ, 728, 85 

Kalas, P., Graham, J. R., Chiang, E., Fitzgerald, M. P., Clampin, 
M., Kite, E. S., Stapelfeldt, K., Marois, C, & Krist, J. 2008, 
Science, 322, 1345 

Krist, J. E., Stapelfeldt, K. R., Golimowski, D. A., Ardila, D. R., 
Clampin, M., Martel, A. R., Ford, H. C, Illingworth, G. D., & 
Hartig, G. F. 2005, AJ, 130, 2778 

Lafreniere, D., Marois, C, Doyon, R., Nadeau, D., Sz Artigau, E. 
2007, ApJ, 660, 770 

Racine, R. , Walker, G. A. H. , Nadeau, D. , Doyon, R. , 
&Marois, C. 1999, PASP, 111, 587 



Leconte, J., Soummer, R., Oppenheimer, B. R., Hinkley, S., 
Brenner, D., Sivaramakrishnan, A., Kuhn, J., Perrin, M. D., 
Roberts, Jr., L. C, Simon, M., Brown, R. A., Chabrier, G., Sz 
Baraffe, I. 2007, in In the Spirit of Bernard Lyot: The Direct 
Detection of Planets and Circumstellar Disks in the 21st 
Century 

Mamajek, E. E., Kenworthy, M. A., Hinz, P. M., & Meyer, M. R. 

2010, AJ, 139, 919 
Marois, C, Lafreniere, D., Doyon, R., Macintosh, B., Sz Nadeau, 

D. 2006, ApJ, 641, 556 

Marois, C, Macintosh, B., Barman, T., Zuckerman, B., Song, I., 
Patience, J., Lafreniere, D., Sz Doyon, R. 2008, Science, 322, 
1348 

Marois, C, Macintosh, B., Sz Veran, J. 2010, in Presented at the 
Society of Photo-Optical Instrumentation Engineers (SPIE) 
Conference, Vol. 7736, Society of Photo-Optical 
Instrumentation Engineers (SPIE) Conference Series 

Marois, C, Racine, R., Doyon, R., Lafreniere, D., Sz Nadeau, D. 
2004, ApJ, 615, L61 

McElwain, M. W., Metchev, S. A., Larkin, J. E., Barczys, M., 
Iserlohe, C, Krabbe, A., Quirrenbach, A., Weiss, J., Sz Wright, 
S. A. 2007, ApJ, 656, 505 

Mesa, D., Gratton, R., Berton, A., Antichi, J., Verinaud, C, 
Boccaletti, A., Kasper, M., Claudi, R. U., Desidera, S., Giro, 

E. , Beuzit, J., Dohlen, K., Feldt, M., Mouillet, D., Chauvin, G., 
& Vigan, A. 2011, ArXiv e-prints 

Metchev, S., Marois, C, & Zuckerman, B. 2009, ApJ, 705, L204 

Mugnier, L. M., Cornia, A., Sauvage, J., Rousset, G., Fusco, T., 

Sz Vedrenne, N. 2009, Journal of the Optical Society of 

America A, 26, 1326 
Oppenheimer, B. R., Kulkarni, S. R., Matthews, K., Sz Nakajima, 

T. 1995, Science, 270, 1478 
Patience, J., King, R. R., de Rosa, R. J., Sz Marois, C. 2010, 

A&A, 517, A76+ 
Perrin, M. D., Sivaramakrishnan, A., Makidon, R. B., 

Oppenheimer, B. R., Sz Graham, J. R. 2003, ApJ, 596, 702 
Pueyo, L., Kay, J., Kasdin, N. J., Groff, T., McElwain, M., 

Give'on, A., Sz Belikov, R. 2009, Appl. Opt., 48, 6296 
Sauvage, J.-F., Fusco, T., Rousset, G., Sz Petit, C. 2007, Journal 

of the Optical Society of America A, 24, 2334 
Serabyn, E., Mawet, D., Sz Burruss, R. 2010, Nature, 464, 1018 
Soummer, R. 2005, ApJ, 618, L161 
Soummer, R., Sz Ferrari, A. 2007, ApJ, 663, L49 
Soummer, R., Ferrari, A., Aime, C, Sz Jolissaint, L. 2007a, ApJ, 

669, 642 

Soummer, R., Pueyo, L., Sivaramakrishnan, A., Sz Vanderbei, 

R. J. 2007b, ArXiv e-prints, 711 
Sparks, W. B., & Ford, H. C. 2002, ApJ, 578, 543 
Thalmann, C, Carson, J., Janson, M., Goto, M., McElwain, M., 

Egner, S., Feldt, M., Hashimoto, J., Hayano, Y., Henning, T., 

Hodapp, K. W., Kandori, R., Klahr, H., Kudo, T., Kusakabe, 

N., Mordasini, C, Morino, J., Suto, H., Suzuki, R., Sz Tamura, 

M. 2009, ApJ, 707, L123 
Wallace, J. K., Burruss, R., Pueyo, L., Soummer, R., Shelton, C, 

Bartos, R., Fregoso, F., Nemati, B., Best, P., Sz Angione, J. 

2009, in Proceedings of the SPIE conference 7440, ed. S. B. 

Shaklan, Vol. 7440 (SPIE) 
Zimmerman, N., Brenner, D., Oppenheimer, B. R., Parry, I. R., 

Hinkley, S., Hunt, S., & Roberts, R., 2011, PASP, 123, 904 
Zimmerman, N., Oppenheimer, B. R., Hinkley, S., Brenner, D., 

Parry, I. R., Sivaramakrishnan, A., Hillenbrand, L., Beichman, 

C, Crepp, J. R., Vasisht, G., Roberts, L. C, Burruss, R., King, 

D. L., Soummer, R., Dekany, R., Shao, M., Bouchez, A., 
Roberts, J. E., & Hunt, S. 2010, ApJ, 709, 733 

Vigan, A. and Langlois, M. and Moutou, C. and Dohlen, K. 2008, 

ASzA, 489,1345-1354 
Thatte, N. and Abuter, R. and Tecza, M. and Nielsen, E. L. and 

Clarke, F. J. and Close, L. M.,2007, MNRAS, 378,1229-1236 



20 



Pueyo et.al 




Fig. 5. — Illustration of the frame classification method. The data consists in a collection of cubes, represented in the middle panel. 
Each cube corresponds to one exposure and each slice within a cube is an image at a given wavelength: here, the temporal dimension is 
represented horizontally and the wavelength vertically. We co-add the frames at the wavelength X po into a target frame T, top panel. We 
create a library of references frames that will be used to build the composite calibrator {Rk}- These references are chosen so that the flux 
of a potential companion at neighboring wavelengths is minimal at the radial location of interest (e.g. minimize self-subtraction). The 
reference ensemble presents both chromatic and temporal diversity. 
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Fig. 6. — Geometry associated with a polar implementation of LOCI. The reference coefficients are chosen in order to obtain the best 
least-squares fit of the target speckle pattern in the large optimization zone, O. The corresponding composite reference PSF is then 
subtracted from the target frame only in a smaller area, the subtraction zone <S, that is included in O. A zoom on the S after LOCI 
subtraction is displayed in the rightmost panel of the figure. 
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Raw Data, H band 



Alcor 



Reduced , J band 



Reduced , H band 



Fig. 7. — Data set used for SNR calculations. Left Column: Raw (top) and processed (bottom) images in the J band. Right Column: 
Raw (top) and processed (bottom) images in the H band. We used a spiral of 11 synthetic companions with a flat spectrum whose angular 
separation vary from 0.3 to 0.5 arcsec. Because of the better performance of the Adaptive Optics system at longer wavelengths, the raw 
contrast is more favorable in the H-band than in J-band. In this example the exposure time is 140 s. Since the average contrast in J-band 
is ~ 10 -4 and ~ 10 -3 in H-band the synthetic companions are detected with a higher signal to noise in the latter case. This figure was 
generated using images reduced by d-LOCI. 
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SNR as a function of NA 
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Fig. 8. — Detectability of faint point sources in the case of an over-constrained LOCI reduction: SNR as a function of the area of the 
optimization zone for constant exposure time and number of references (here the exposure time is 140 s and N-jz < 23). The synthetic 
companions illustrated on Fig.PFlwere used for these calculations: the figure shows the average of the post LOCI SNR over the 11 companions. 
Top curves: H-band results. This shows that in the over-constrained regime d-LOCI performs almost as efficiently as classical LOCI as 
far as the detectability of faint point sources is concerned. While these results were obtained using g = 1, reduction with g ^> 1 yield the 
same qualitative behavior, albeit at slightly higher SNR levels. Bottom curves: J-band results. Details about the second order trends in 
this figure are discussed in the body of the present manuscript. 
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SNR as a function of number of cubes 
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Fig. 9. — Detectability of faint point sources in the case of an under-constrained LOCI reduction: SNR as a function of the number 
of cubes reduced (increasing exposure time and increasing number of references), for a constant area of the optimization zone, constant 
Na- The synthetic companions illustrated on Fig. [7] were used for these calculations: the figure shows the average of the post LOCI SNR 
over the 11 companions. In the case of under-constrained problems, the relative performances of classical LOCI and d-LOCI with respect 
to detectability are a function of the contrast of the companion. Top curves: H-band results. For bright companions, H-band results, 
classical LOCI provides better noise suppression. Bottom curves: J-band results. For fainter companions, J-band, the noise suppression 
with classical LOCI is associated with a large signal depletion (e.g, the companion's signal is fitted by one of the references). When using 
d-LOCI, this depletion is mitigated, and thus larger SNR are obtained. While these results were obtained using g = 1, reduction with 
g ^> 1 yield the same qualitative behavior, albeit at slightly higher SNR levels. 



